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Abstract: Many multi-loop calculations make use of integration by parts relations to 
reduce the large number of complicated Feynman integrals that arise in such calculations 
to a simpler basis of master integrals. Recently, Gluza, Kajda, and Kosower argued that 
the reduction to master integrals is complicated by the presence of integrals with doubled 
propagator denominators in the integration by parts relations and they introduced a novel 
reduction procedure which eliminates all such integrals from the start. Their approach has 
the advantage that it automatically produces integral bases which mesh well with gener- 
alized unitarity. The heart of their procedure is an algorithm which utilizes the weighty 
machinery of computational commutative algebra to produce complete sets of unitarity- 
compatible integration by parts relations. In this paper, we propose a conceptually simpler 
algorithm for the generation of complete sets of unitarity-compatible integration by parts 
relations based on recent results in the mathematical literature. A striking feature of our 
algorithm is that it can be described entirely in terms of straightforward linear algebra. 



Keywords: [NLO Computations . 



Submitted to JHEP 



Contents 



[l]. Introduction 



|2|. Preliminaries 



||. The Algorith m 







An Explicit Example 



11 



|5]. Conclusions 



17 



1. Introduction 

When the technique of integration by parts in d dimensions was first proposed by Tkachov 
and Chetyrkin Q, it represented a major breakthrough in the study of perturbative gauge 
theories at the multi-loop level. 1 Their discovery was of fundamental practical importance, 
as it allowed researchers to perform many multi-loop calculations that were previously 
thought to be intractable. Furthermore, the idea of the technique is simple to describe. 
By taking integrals of derivatives in d dimensions, one generates a tower of equations for 
the Feynman integrals belonging to a particular topology. Then one tries to solve these 
equations, either by inspection or via some systematic procedure, for an independent basis 
of master integrals. 

Unfortunately, the total number of equations generated in this way grows rapidly with 
the number of loops and external states in the integral topology under consideration. As 
a consequence, the solution of the so-called integration by parts relations is complicated 
for all but the simplest examples. For many years after the technique was introduced, no 
systematic procedure for the solution of integration by parts relations was known and it was 
therefore only possible to apply the method to simple integral topologies. 2 The situation 
changed just over a decade ago with the introduction of a Gaussian elimination-like solution 
algorithm. This algorithm, due to Laporta [||], was a crucial step forward because it allowed 
researchers to apply the integration by parts technique to highly non-trivial problems for 
which an ad hoc approach would be impractical if not impossible. Although Laporta's 
algorithm has been tremendously successful, it has long been known that it may lead to 
master integrals with doubled propagator denominators. 

To be more precise, for an L-loop topology, let V = {£i, ■ ■ ■ ,£l, ki . . . , &jv} be the set 
of loop momenta together with a set of N independent external momenta and let v be a 
generic element of this set. Normally, one generates the set of integration by parts relations 
in an obvious way, considering each v 6 V in turn and writing all possible equations of the 

1 For the reader less familiar with integration by parts, we strongly recommend Smirnov's excellent book 
on Feynman integral calculus B . He explains the technique and works through a number of simple (and 
non-simple) examples. 

2 Of course, simple is a relative term. For a rather impressive early application of the method at the 
three- loop level see the long write-up for MINCER H. 
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form 

n = T\( [ JlIl) JL p' ) ^i,-4) ai -^i,..4) tt A , . 

11 \J (2*)*) dij'\ D 1 (£ 1 ,...l L fi---D m (£ 1 ,...£ L )b™ )■ [ ■ > 

Crucially, the indices 04 and bi satisfy certain boundary conditions; typically the irreducible 
numerators (here we assume that each Mi has mass dimension two) for a given topology 
enter raised to, at most, some relatively small positive integer and, in some cases, certain 
propagator denominators are constrained to have non-negative indices since otherwise the 
resulting integrals vanish in dimensional regularization. This set of equations together with 
the appropriate boundary conditions can then be fed into Laporta's algorithm. The price 
one pays for being able to trivially generate the set of integration by parts relations in this 
fashion is twofold. Not only will the algorithm (Laporta or typical variation thereof []|, || [7j) 
used to solve the resulting system of equations typically have to eliminate an enormous 
number of spurious integrals as it attempts to solve the system, it is not straightforward to 
ensure that each integral in the basis ultimately returned by the algorithm has the property 
that b{ = 1. 

Although there is nothing wrong with integrals that have some bi > 1, they may 
be inconvenient for particular multi-loop applications. For example, in computational 
approaches based on generalized unitarity, one would like to have a basis with well-defined 
unitarity cuts in all channels (see e.g. || and |J for recent work in this direction at two 
loops). It is unclear how one would make sense of integrals with doubled propagator 
denominators in such a framework. Gluza, Kajda, and Kosower (hereafter referred to as 



GKK) also argued in [1C] that master integrals with doubled propagator denominators can 
be significantly harder to expand in e than those without. 

The idea of the GKK procedure is relatively easy to state now that the stage is set. 
With the above motivation, GKK found that they could completely avoid the introduction 
of doubled propagators by imposing bi = 1 from the beginning. They observed that, 
generically, one expects doubled propagators for the simple reason that the derivatives in 
eqs. ( [l . l| ) act on the propagator denominators, D^. They also recognized that there is no 
good reason why one ought to consider the elements of V one at a time; one can generalize 
eqs. (|l . l| ) by replacing u^) with a general linear combination of the elements of V: 

yt([£u\ g ( E^i N ^ ) vi j) ^i(h,...£Lr---M q (£ 1 ,...£ L rA 

f = \\J (27r)«V d£, { D 1 (£ 1 ,...£ L fi...D m (£ 1 ,...£ L fm J ' ^ ' } 

It is also convenient at this point to combine together some of the equations by summing 
over j: 




fi" a?'^? Mi(h, . . . £ L ) ai ■ ■ ■ M q (h, . . . i L r« ' 



; , , v . , ..£ L )^...D m (h,...£ L f 

In a nutshell, the GKK strategy is to start with all the bi equal to one 

L ' r d d £i \ A d (Y:^ afv\ i] Nx{£x,...£L) a ^--M q {h,...£ L ) a « 



■ (1.3) 



11 V (2vr)V j^dlj \ D 1 (£ 1 ,...£ L )---D m (£ 1 ,...£ L ) J { - ) 
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and then choose the coefficients a\ in such a way that the numerator exactly cancels all 
unwanted, derivative-generated powers of the propagator denominators. In other words, 
for each k, we demand that 

L L+N 

E E a i ] -Qf ■ v ? ■ ■ ■ i L ) ai ■ -Wi, ■ ■ ■ i L ) a " = -2T^-^D k , (1.5) 

j=i i=\ i 

where T^ k,aiy " ,a< ^ is some polynomial built out of the irreducible numerators and Lorentz 
invariant combinations of the vectors in V? Each independent set of L(L + N) coefficients 
&i satisfying eqs. (|l.5D , upon substitution into eqs. fll.4|) , yields a unitarity-compatible 
integration by parts relation, by construction free of doubled propagators. 

The downside of this novel approach is that now one needs to find a way to generate 
complete sets of a^p coefficients and it turns out that this is not a straightforward task. 
Although GKK did propose a solution 4 to this problem in p0|| , they admitted that their 
solution was somewhat provisional and that there is likely room for improvement. The 
GKK solution relies heavily on the use of Grobner bases, important constructs in com- 
putational commutative algebra which are, however, notoriously difficult to compute in 
practice M]. 5 In this paper we propose a completely different solution to the problem 
of generating complete sets of unitarity-compatible integration by parts relations. As we 
shall see, our Algorithm [l] is based entirely on simple linear algebra and, in particular, 
completely avoids the use of Grobner bases. 

This article is organized as follows. In section ^ we describe more precisely the problem 
we wish to solve and introduce some useful notation. In section ^ we present the main result 
of this paper, Algorithm [l], and talk the reader through it. In Section |], we give a detailed 
example of how our algorithm works in practice. In section [5] we present our conclusions 
and outline our plan for future research. 



2. Preliminaries 



In this section we take a closer look at the reduction procedure introduced by GKK and 
discuss its strengths and weaknesses. Our initial goal will be to precisely set up the math- 
ematical problem that lies at the heart of the GKK procedure and discuss why (in the 



opinion of both GKK and the present author) the solution presented in [10] is not likely to 
be the best one possible. We will then explain our approach to the problem and illustrate 
with a very simple example what precisely Algorithm [l] is designed to do. 



3 We write _2r (fc ' ai ' , ,,a<j) instead of r( fc >°i.— ,°- q ) so t ^ at we w jrj ultimately arrive at the same form that 
GKK did in their eqs. (4.1). 

4 Actually, GKK presented two distinct algorithms for the generation of complete sets of unitarity- 
compatible integration by parts relations. In what follows, when we refer to GKK's "solution" or "solution 
algorithm" with no additional qualifier, we are referring to their best solution (what they call Algorithm 

in). 

For the reader less well- versed in computational commutative algebra, we can strongly recommend the 
very well-written and concise survey by Adams and Loustaunau 113] . Most of the relevant mathematical 
concepts are also defined and briefly explained in the GKK paper [tl0|. 



-3- 



As we left them, eqs. ( |1.5| ) look rather cumbersome. We can clean them up by 
absorbing the numerator polynomial N\{1\, ■ ■ ■ (-lY 1 • • • M q {i\, ■ ■ ■ ^l) " 1 into the alp 

L l+n „_ 

j=i i=i 3 

with the understanding that now the ccp are dimensionful. We will have to take this into 
account in our search for independent solutions. To make further progress, GKK expressed 
( |2.1| ) as a matrix equation: 

P-E = 0, (2.2) 

where 

= (a[ l) , . . . , a^ +N , , a[ L) , . . . , a^ N ,T^ , . . . , T^) . (2.3) 

Given an explicit expression for V and an ordering on the set of propagator denominators, 
the entries of E can be straightforwardly read off from (|2.l| ). For an explicit example, we 
refer the reader to eq. (5.3) of |I(J. In eq. (5.3), GKK wrote out E explicitly for the planar 
massless double box, choosing - 2l\ ■ k\,t\ - 2l\ ■ k\ - 2£± ■ fc 2 + si2,^i,^| - 2^ 2 ■ 

&4, £2 — 2^2 • ks — 2^2 • ^4 + S12, i\ + £\ + 1l\ • h} for the ordering on the set of propagator 

denominators (this ordering fixes the sequence of the columns of E). 6 The point of making 

(i) 

this rearrangement is that now the problem of determining all independent sets of a\ 
coefficients looks like a well-known, well-studied problem from computational commutative 
algebra: the computation of the syzygy module 7 of a module for which one has an explicit 
set of generators. 

Actually, as explained by GKK, it suffices to solve this problem for ideals, since given 
a set of generators for a module A4, one can easily define a set of generators for an ideal 
with exactly the same syzygies. Suppose the set {Mi, . . . ,M r } generates M.. By simply 
taking the dot product of each Mj with a tuple of dummy variables satisfying the relations 
Utj = 0, {t±, . . . ,t m }, we can convert our generating set for M. into a generating set for 
an ideal X canonically associated to Ad. The class of ideals canonically associated to the 



modules described by eqs. (2.1) can be taken to be homogeneous of uniform degree two. 
We can see this as follows. By definition, each parameter entering into the matrix E of eq. 
(|2.2| ) is of the form • or and is going to have mass dimension two. Therefore, 
it makes sense to take every independent Lorentz product that can arise in E to be an 
independent variable and the class of modules described by eqs. ( |2.1| ) to be homogeneous 



of uniform degree one. It then follows, after applying the canonical map described above, 



6 We should also point out that the precise definitions GKK made for the of (what they call ) and 
T' fe ' a ) (what they call Uk) do not seem to exactly reproduce eq. (5.3). However, the differences can be 
taken into account by appropriately rescaling the unknowns and are therefore unimportant. 

7 Given a set of generators for a module M, {Mi, . . . , M r }, the syzygy module of M is simply the set 
of all P = (pi, . . . , fe r ) such that ' Mi = 0. In this paper, we will deal only with submodules of 

F[xi, . . . , x n ] r and F[xi, . . . , x n ] m for a field, F. 
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that the ideals of interest are homogeneous of uniform degree two. In the end, we find that 
each generator of T has terms of the schematic form atj x&, a £ F for some field F. 8 

Before continuing, we need to introduce a little more notation and point out an impor- 
tant fact about the syzygies of generating sets for homogeneous ideals of uniform degree. 
Given an ideal generated by P = (pi, . . . ,p r ), we define Sd(P) to be the set of all syzygies 
of P such that /3 = (61, . . . ,b r ) € Sd(P) implies that each 6j is a polynomial of uniform 
degree d. We will typically refer to Sd(P) as the set of all degree d syzygies of P. It turns 
out that, for the case of homogeneous ideals of uniform degree, the syzygy module, S(P), 



is a graded module [|13[: 



S(P) = QS d (P). (2.4) 

d=0 

This means that, for homogeneous ideals of uniform degree, it suffices to search for degree 
d syzygies. 

Determining a basis for the syzygy module of a generic ideal is known to be a very 
difficult problem (and is very much an active area of mathematical research). Therefore, 
one needs a dedicated solution algorithm, tailor-made for the class of ideals described above. 
As mentioned in the introduction, the solution algorithm employed by GKK relies heavily 
on the use of Grobner bases. GKK chose the well-known Buchberger algorithm to 
compute their Grobner bases. In their paper, GKK pointed out that there have been 



a number of recent attempts to improve on Buchberger's algorithm (see e.g. [15] for a 
description of one of the most promising of these recent attempts, based on Faugere's F5 
algorithm and concluded that their Buchberger-based approach was not likely to be 
optimal. In this paper we rethink their approach at a more fundamental level. 

Certainly, one could attempt to compute complete sets of syzygies using an approach 



based on Faugere's F$ algorithm [17] or some other improved algorithm for the computation 
of Grobner bases. However, it is actually no longer clear that one should use Grobner bases 
at all. Quite recently (after the appearance of |i~0| l) it was shown in |l3| that, remarkably, 
bases for the modules of syzygies of special classes of ideals can be computed using simple 
linear algebra. Actually, at this stage, the reader may be wondering what stopped GKK 
from using a linear algebra-based approach in the first place. Naively, linear algebra seems 
to offer a very easy way to compute syzygies. 

To illustrate why the ideas of reference [jDJ are non-trivial, we consider the ideal 
generated by P = {x\ — X2,2x2 — xi} and attempt to compute its syzygies using linear 
algebra. 9 By inspection, we see that P has no degree zero syzygies. By definition, a 
degree one syzygy of P must have the form {c\X\ + 02X2,03X1 + C4X2) for some elements 
(c\, C2, C3, C4) of Q. Starting with this ansatz, we can take the dot product with P 

(ciXi+C 2 X 2 , C3X1+C4X2) P = (ci-C3)Xi + (-Ci+C2 + 2c3-C4)xiX2 + (-C2 + 2c4)x2 , (2.5) 



3 In this paper, F is given by the field of rational functions of the auxiliary parameters in the matrix E 



{e.g. xu in the case of the planar massless double box treated in detail by GKK) of eq. (2.2). However, 
for practical purposes, it is more convenient to simply assign prime numbers to the dimensionless auxiliary 
parameters and work over Q. 

9 For the sake of definiteness, suppose that we are working with polynomials in Q[xi, 052]. 
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set the coefficient of each monomial in the sum to zero, and solve the resulting system 
of equations. It turns out that there is a one-parameter family of solutions which we 
parametrize by C4: 

(C1,C 2 ,C 3 ,C 4 ) = (-C4,2C 4 ,-C4,C 4 ) . (2.6) 

We can arbitrarily set C4 = 1 to obtain a basis for Si(P), (— xi + 2x2, —%i + X2). So far 
so good. By definition, a degree two syzygy of P must have the form {c\x\ + C2X1X2 + 
C3X2, 04x^+05X1X2+06X2). Going through the same procedure, we arrive at a two-parameter 
family of solutions which we parametrize by C5 and c§: 

(C1,C2,C 3 ,C4,C5,C 6 ) = (-c 5 - c 6 ,2c 5 + c 6 ,2c 6 , -c 5 - c 6 ,c 5 ,c 6 ) . (2.7) 

By replacing (05, cq) with each of the standard basis vectors for M 2 in turn, we find that 
(— 1, 2, 0, —1, 1, 0) and (—1, 1, 2, —1, 0, 1) form a basis for the set of solutions. In order to 
map these solutions back to degree two syzygies of P, we partition each vector into two 
non-overlapping subsets of length three without changing the overall ordering of the entries 

{ ({-1, 2, 0}, {-1, 1, 0}) , ({-1, 1, 2}, {-1,0, 1}) } (2.8) 

and then dot each resulting 3-tuple into {x\, x\X2, x|}: 

{x 2 ,xix 2 ,X2} • {-1,2,0},{x 2 ,xix 2 ,X2} • {-1,1,0V 
^{x 2 ,xix 2 ,X2} • {-1,1,2},{x 2 ,xix 2 ,X2} • {-1,0,1})| 
= I ( - x\ + 2xix 2 , -xl + X1X2J , ( - xl + xix 2 + 2x3, ~ x l + x i) } • ( 2 - 9 ) 

This overly formal description of the mapping back to syzygies of P is overkill for the 
example at hand but will be useful later on. 

We might be tempted to conclude that we have found two new, linearly independent, 
degree two syzygies of P. This, however, is not true. It turns out that the syzygies of eq. 
Q2.9[) can be expressed as multiples of (— x% + 2x2, — %i + X2). Explicitly, we have 



{ 



( — x\ + 2xix 2 , -xl + xix 2 ) = xi(-xi + 2x 2 , -xi + x 2 ) ( 2 -10) 

and 

^ — x\ + xix 2 + 2x|, ~ x i + x t) = ( x i + x 2){— x\ + 2x 2 , -xi + x 2 ) • (2-11) 

Besides highlighting the profound difference between linear independence in a vector 
space and linear independence in a module, this example shows what goes wrong if one tries 
to compute the syzygies of homogeneous ideals of uniform degree using naive linear algebra. 
One is able to compute syzygies in a straightforward manner but not a basis of linearly 
independent syzygies. It is therefore remarkable that, under appropriate assumptions, this 
obstacle is actually surmountable. In the next section, drawing upon some of the ideas 



introduced in [ 13 1 , we present a simpler linear algebra-based alternative to the solution 



adopted by GKK. 
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3. The Algorithm 



The purpose of this section is to give an explicit pseudo-code detailing our solution to 
the problem defined in Section |2| and to thoroughly comment it. The pseudo-code we 
present (Algorithm [j] and associated subroutines) is quite explicit and should allow the 
reader to fashion a crude implementation of our algorithm in Maple or Mathematica with 
very little effort (beyond that required to understand the algorithm in the first place). We 
should emphasize that we do not in any way claim that an implementation based on our 
pseudo-code is an effective implementation. The pseudo-code given below is intended to 
be maximally clear as opposed to maximally efficient. 10 Before discussing the non-trivial 
features of Algorithm |l[ we need to make a few more definitions. 

First, let M^ 1 ' = Tx[ , . . . , X ^j_ n _ 1) , j be the set of all monomials of degree d built 

d!(n-l)! - 1 

out of the n variables {xi, . . . , x n }. u For example, = {x\ , X1X2, xfy- For definiteness, 
we will always order sets of monomials lexicographically. This choice of ordering is just a 
choice and has no deeper significance. Now, if we have in hand a sequence of r homogeneous 
polynomials of uniform degree two, Pq = {p\, . . . ,p r }, then Pd is defined to be the outer 
product of and Po : 

P d = { Pl X[ d) piXg) , , Pr X{ d) p r X {d l n _ iy , }. (3.1) 

k d!(n-l)! d!(n-l)! J 

Clearly, Pd is a set of r ^f^" nr homogeneous polynomials of uniform degree d + 2. For 
example, if we take Pq = P = {xi— X2, 2x2— %i} it is easy to see that P2 = {(x\— X2)x\ , (x\ — 
X2)x\X2, {x\ — X2)x\^ (2x2 — xi)xf, (2x2 — x\)x\X2 , (2x2 — x\ )x|} . The construction of P2 can 
be thought of as an intermediate step towards the extraction of the degree two syzygies of P. 
Instead of solving the system (01x^+02X1X2+03X2, c^x\ +05X1X2+06X2)- {xi— X2, 2x2— xi} = 
0, we can construct P2 and then solve the system (ci, 02, 03, 04, 05, cq) • {(xi — X2)xf , (xi — 
X2)xiX2, (xi — X2)x|, (2x2 — x\)x\, (2x2 — xi)x\X2, (2x2 — ^1)^2} = 0- So, instead of trying 
to work with Sd(Po) directly, we can work with the vector space So(Pd). This natural 
correspondence between basis vectors of So(Pd) and degree d syzygies of Po is an essential 
part of Algorithm [j], which is why we took the time to carefully describe it while working 
through the illustrative example at the end of Section ||. However, it is important to 
remember that this map does not necessarily yield independent elements of S(Pq); to work 
as advertised, our solution algorithm must be able to determine, using only linear algebra, 
what basis vectors of So(Pd) correspond to new syzygies of Poj linearly independent of 
those already determined. 



For example, the first statement in the upper branch of the If statement in Algorithm |l| is completely 
superfluous and was included only because, in our opinion, it makes the functionality of Subroutine |l| 
significantly easier to understand. 

11 Note that there are precisely ^fe^ monomials of degree d built out of n variables. This result follows 
immediately once one recognizes that the monomial counting problem is isomorphic to one of the usual 
"ball-box" counting problems of enumerative combinatorics. We refer the unfamiliar reader to Section 1.4 
of Stanley's textbook on the subject ]lq] . 
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Algorithm 1 

Require: A minimal generating set for a homogeneous ideal of uniform degree two, Pq = 
{pi(h, ...,t 

mj . . . , x n ), • • • , Pr(ti, ■ ■ ■ , t m , X\, . . . , x n )}. The pi are further constrained to 
have terms of the schematic form atjXk with a £ F and ttfj = for all i and j. 
SyzListg = Bq = {0}. For a given integral topology with generic irreducible numerator 
Afi(£i,... i L ) ai ■ ■ ■ Af q (£i, . . . Il)"" 1 , A is the largest value of Ya=i a i allowed by the bound- 
ary conditions on the Oj (we assume that each Mi has mass dimension two). M^ 1 = 
<x[ d \ . . . ,^[d+n-i)! } is the lexicographically ordered set of all monomials of degree d built 

d\(n-l)\ > 

out of the n variables {x%, . . . , x n }. 



•Set d = 1 
while d < A do 
if B d _ x ^ {0} then 

(n) 

•Take the outer product of Po and M^J-, to determine P d -i' 

p _ /_ _ y(<i-l) ' _ v (*-l) _ \ 

Pd-l — \Pl A l , ...,PlA (d+w _ 2)! , ,---,Pr^ (d+ n -2)> ( 

(<j-l)!(n-l)l (d-l)!(n-l)! J 

(n) 

• Take the outer product of Po an d AO to determine P d : 

Pl _L _ y W _ _ \ 

, ■ ■ ■ ,PlA (d+n _ 1) , , ,Pr^\ r--iM (j+n-l)! f 

d!(n-l)! d!(n-l)! J 

•Apply Subroutine [l] to each element of P^-i and call the union of the results A: 

\B d -i\ 

A= (J <7(ai) 

8=1 

•Regard the elements of A as the rows of a matrix, find a row echelon form of this 

matrix, keep only rows with at least one non-zero entry, and call the result C 
•Replace with zero the s entries of Pd which correspond 

to the pivot columns of C and call the result G 
•Apply Subroutine ^ to G and call the result D: 
D = SSyz(G) 

•If D ^ {0} then B d = C |J D and if D = {0} then B d = C 
•Apply Subroutine || to D to determine SyzListg : 
SyzListg = ToSyz(D) 

else 

In) 

• Take the outer product of Po and M d to determine P d : 

p _ /_ y(<9 _ yM _ _ \ 

, ...,piA (d+ „_ 1)! , ,Pr*i ,PrA (d+n-iy. f 

k d!(n-l)! d!(n-l)! J 

•Apply Subroutine || to P^ and call the result D: 

D = SSyz(P d ) 
•B d = D 

•Apply Subroutine || to D to determine SyzList d : 
SyzListg = ToSyz(D) 
end if 
•d = d + l 
end while 

return f [j SyzList d ^ ^{0} 



-8- 



Subroutine 1 a (a) maps a G £?<z_i to n linearly independent elements of So(Pd) 



Require: The objects introduced in Algorithm |l] and a = yA\,...,A r 



r(d+n-2)\ 



(d-l)'.(n-iy. 

A m xr^ + ...+A_ r(d+n-2)\ P r X (d+n _ 2)! G Bd_l 

(d-l)!(n-l)! (d-l)t(n-l)! 



return <MipiX^ ^xi H \- A r (d+n-2y. p r X^ d (d ]} n _ w x±, ,Aij>iX{ d X 'x n + 

(d-l)!(n-l)! 

Pr^ ( d+n _2)i x n f 

(d-l)!(n-l)! (d-i)!(n-l)! J 



(d-l)!(n-l)! (d-l)!(n-l)! 

/j 1 

r(d+re-2)! 



Subroutine 2 SSyz(G) computes a basis for Sq(G) 

Require: The objects introduced in Algorithm [l] and a sequence, G = {gi, . . . , q n }, of homoge- 
neous polynomials in the variables {t%, . . . ,t m , x%, . . . , x n } with coefficients in F. 

• Introduce a sequence of |G| dummy variables: 

y = {vi,'-,v\g\} 

•Expand the dot product G ■ Y and partition the terms in the sum into equivalence classes, 
declaring terms equivalent if and only if they contain the same monomial in the variables 
{ii , . . . , t m , x\ j . . . , Xfi\ 

•Set the sum of terms in each equivalence class equal to zero, collect the resulting set of 
equations constraining the y{, and drop any equations that are obviously redundant 

•Put the remaining equations into matrix form: 
QY = 

return a basis for Null(Q) 



Subroutine 3 ToSyz(D) maps v elements of So(Pd) to v elements of Sd(Po) 

Require: The objects introduced in Algorithm |] and a list, D = {d\, . . . , d u }, d{ G So(Pd)- 

•Partition each vector d{ into r distinct non-overlapping sequences of length ^z^rw > 
taking care to preserve the order of the entries of each df 
*-►{*?>,. ..,4°} 
return { (ilff - £\ . . . , ■ , , (m™ ■ . . , ^ • *<">) } 

We now describe how Algorithm [l] works in some detail. By assumption, the m prop- 
agator denominators of the integral topology under consideration form a linearly indepen- 
dent set. 12 This implies that Pq is a minimal generating set for the ideal under consideration 



12 If this is not the case, one should first reduce to simpler integral topologies satisfying this condition 
before attempting to apply Algorithm [[]. GKK explain this procedure in some detail in Jl(i| |. 
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and there are no non-trivial syzygies of degree zero. Since Algorithm [l] proceeds incremen- 
tally in the degree of the syzygies, it is convenient to define two bookkeeping lists, and 
SyzList rf , indexed by d. Bj is a basis for the vector space So (Pi) and SyzList^ is a set of 
linearly independent elements of S^(Po) which are, for d > 1, also linearly independent of 
all syzygies of degree < d — 1 in the set SyzListj (we will call the elements of SyzList rf 
new degree d syzygies of Po). Due to the fact that Pq has no degree zero syzygies, both lists 
are initialized to {0} and d is initialized to one. Since Bo = {0}, Algorithm [l| always starts 
with a pass through the lower branch of the If statement. Using Subroutine Algorithm 
|l] determines a basis for So (Pi) along the lines described in the example near the end of 
Section ||. For the sake of discussion, we assume that this basis is non-trivial. This basis, 
D, is then used to determine both B\ and, via Subroutine ||, SyzListj. In this section we 
will not step into Subroutines || or |3| There is nothing non-trivial about them and we will 
in any case go through the subroutines once explicitly in Section f|. 

Increment d to two. Subroutine || lifts the syzygies in B\ to syzygies in So(F^) in every 
independent way possible. This is accomplished by simply multiplying each of the elements 
of B\ (written out explicitly in terms of their coordinates, the elements of Pi) by each one 
of the variables {x\, . . . ,x n } in turn. The results can then be interpreted as syzygies in 
5*0 (P2) (written out explicitly in terms of their coordinates, the elements of P2). Clearly, 
if /3 is a scalar syzygy of Pi , f3 ■ (xj Pi ) =0 for each 1 < i < n by linearity. This implies 
that /3 can be interpreted as a scalar syzygy of P2 in n independent ways. Given v syzygies 
in So (Pi), Subroutine [l] produces vn syzygies in So(P2) (some of which may be linearly 
dependent). As we shall see, Subroutine [l] is a crucial first step towards determining which 
syzygies in So(P2) correspond to new degree two syzygies of Po and which do not. 13 

Since the syzygies in So(P2) produced by Subroutine |l| are not guaranteed to be linearly 
independent, the next step is to put the output of Subroutine [l] into row echelon form, 
discard all rows without non-zero entries, and call the result C . Since, by construction, 
each of the s rows of C is a scalar syzygy of P2, we could in principle rewrite s of the 
polynomials in P2 as linear combinations of the other r ^j^_~^ ! — s. Algorithm [l] takes this 
fact into account by replacing the s entries of P2 which correspond to the pivot columns 
of C with zero and calling the result G. In constructing G, it has isolated the subspace 
of So(P2) which is in correspondence with the new degree two syzygies of Po (if any new 
degree two syzygies exist). Next, Algorithm |l| uses Subroutine [2] to determine a basis, 
D, for the subspace of So(P2) under consideration. However, in order to search for degree 
three syzygies in an analogous fashion, we need a basis for So(P2) in its entirety. Therefore, 
Algorithm |l| sets B2 = C [J D (for the sake of discussion assuming that Subroutine ||| found 
new degree two syzygies). Finally, Algorithm [l] solves for the new degree two syzygies 
themselves by applying Subroutine || to D. After incrementing d to three, Algorithm || 
would pass through the upper branch of the If statement again in an attempt to find new 
degree three syzygies of Po. 

13 The reader who has read reference [jL3j might be concerned that we have not yet spoken about the 
principal syzygies of our ideals. Actually, the fact that the ti are simply coordinate variables with no 
independent existence of their own (they satisfy titj = for all i and j) implies that the ideals of interest 
to us have no principal syzygies at all. 



- 10 - 



Of course, we have been assuming throughout our discussion that A > 3. It is worth 
emphasizing that, the termination condition we are using for Algorithm || comes entirely 
from the physics. The fact that irreducible numerators typically have, at worst, a relatively 
small mass dimension is the only reason that we were able to fruitfully apply the ideas of 
reference [13| and construct Algorithm |l]. Although we believe that the treatment given 
here is appropriate, some readers may prefer a more formal one. If that is the case, then 
we recommend reading [13]. Most of the non-trivial aspects of our pseudo-code are treated 
there as well in the style preferred by mathematicians. 



4. An Explicit Example 

In this section, we show how Algorithm |l| works in practice by going through it for a 
particular example. It was challenging to find a physically motivated example compact 
enough to present in detail and, at the same time, rich enough to give the reader a good 
feeling for how the algorithm functions. In the end, we found that the module given by 
the irreducible part of the module associated (associated in the sense of the construction 
reviewed in Section ||) to the planar massless double box works very well. In the solution 
algorithm of GKK, the study of this module (hereafter referred to as A4) is the first step 
towards the determination of the complete set 14 of linearly independent syzygies of the 
module associated to the planar massless double box. In their paper, GKK assert that A4 
has just three linearly independent syzygies: one of degree one and two of degree two. This 
a priori knowledge of the syzygy module of A4 will allow us to stop the example when 
it ceases to be interesting. Otherwise we would have to make several more trips through 



the While loop (as explained in [1C], A = 6 for the planar massless double box), each 
time discovering no new syzygies. For the sake of clarity, our exposition will mirror the 
pseudo-code of Section |3| quite closely. 

As a preliminary step we must derive the generators of A4 . In Section |2], we pointed 
out that the ordering adopted by GKK for the propagator denominators of the planar 
massless double box is given by {t\,l\ — 2£\ ■ k\,£\ — 2£\ ■ k\ — 2£\ ■ hi + sviA%% ~ 
2£ 2 ■ k 4 Jl - 2£ 2 ■ h - 2£ 2 ■ k 4 + s 12 ,l\ + l\ + 24 • £ 2 }. We actually prefer the ordering 
{£\ , £\ - 2h ■k 1 ,£ 2 1 -2£ v k 1 -2£ v k 2 +s 12 ,£ 2 2 - 2£ 2 ■ k 3 -2£ 2 ■ k 4 +s 12 , l\ - 2£ 2 -k^£\,£\ +£ 2 2 +2£ x -£ 2 } 
and this is what we will use. Note, however, that we do adopt their ordering for the 
generators themselves. By definition, the generators of A4 are the generators of the module 
associated to the planar massless double box by eqs. (pll) (the rows of the matrix E in eq. 



Here it is perhaps worth pointing out that, in fact, the application of our algorithm to the module 
associated to the planar massless double box yields more independent solutions than GKK found (working 
modulo reducibility as they do). We conjecture that, perhaps, GKK did not really seek the complete set of 
linearly independent syzygies modulo reducibility but were instead content to determine a subset sufficient 
for the elimination of as many irreducible numerators as possible. If our reading of GKK is correct, then 
the obvious question is whether discarding potentially useful linear relations between Feynman integrals of 
a given topology is prudent. We suspect that this is not the best strategy because, at least for the planar 
double box, the elimination of irreducible numerators in this manner seems to lead to a large number of 
irreducible integrals of simpler topology. Unfortunately, further exploration of this interesting question is 
beyond the scope of the present paper. 
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(5.3) of |10|| ) reduced over the propagator denominators of the massless double box. This 
reduction is effected by making the substitutions {/C3 — > — k\ — k 2 — k^, if — > 0,4 ■ k\ — > 
?i - 4 ^ 0,4 ■ k 2 

A4 is generated by 



0,£l ->0,£ 2 -h-> 0,4 -e 2 ->0,£i-k2-> s 12 /2,£ 2 ■ k 2 -> -£ 2 ■ ki - s 12 /2}. We find that 



0,0,-^,0,0,0,0 



( S\ 2 

V 2 • 



0,0,0,0,0,-4 ■ h 
0,0,0,^,-4^4,0,0 

+4-4,4-4,4-4), 

0,0,0,-^,0,0,4 -k 4 



0,-4-4,^,0,0,0,0 



0,0,-^,0,0,0,4 -h 



^.^,-^^ + 4.^,^+4.^,0,0,0,4^4) . 



(0,0,0,-^,0,0,0 



0,0,0,-4 • h 



X14S12 



n n n Sl2 1 1 X14S12 

0,0,0, — + 4 • 4, - 

4 • ki, - 4 ■ h, -h ■ h 



(4.1) 



If we let x% = 4 • ki, x 2 = 4 • 4, £3 = S12, and take the dot product of each generator 
in ( |4.1| ) with (4,4,4,4,4,4,4), we arrive at the generating set for the ideal (hereafter 
referred to as I) canonically associated to A4, Pq: 



Pa 



1 4^3 4^3 4^3 

2*3^3, -5 t 2 X 2 , t 7 X 2 — , — 4x 2 , 



1 



:4^3Xi4 + hxi + t 2 X! + 4x1 



4^3 4 X 3 l 

+4^1 + -z- , — ; t 5 Xi, --4^3, 



1 4^3 
'2*5^3X14 + UX 2 + 4^2 + UX 2 + t 7 X 2 + 



1 



:4a;3X:i4 - 4^2 - 4^2 - 4^2 - 4^2 



4^3 



,4x1 



4x3 1 



(4.2) 



Only the first ten rows of E remain non-zero after the reduction is carried out. Provided 
that all m propagator denominators are independent of one another (which is certainly 
true in our case), the Po we arrive at in this fashion will always be a minimal generating 
set for I which implies that, as assumed in Algorithm [l], SyzListg = Bq = {0}. Actually, 
for most of steps of the algorithm the explicit form of Pq is unimportant and we suppress 
it, writing instead P = {pi,p 2 ,P3,P4,P5,P6,P7,Ps,P9,Pio}- 

We initialize d to one and enter the While loop. Since Bq = {0}, the algorithm directs 

(3) 

us to the lower branch of the If statement. M x = {x\, x 2 , X3}, Pq = {pi,p 2 ,Ps,P4,P5,Pe,P7, 
P8,P9,Pio}, and 



Pi = {PIXI , p\X 2 , piX 3 , p 2 X\ , p 2 X 2 , p 2 X3 , p 3 Xl , p 3 X 2 , P3X3 , P4X\ , P4,X 2 , P4X3 , p 5 Xi , p 5 X 2 , P5X3 , 
p 6 Xi , p 6 X 2 , P6X3 , P7X1 , p 7 X 2 ,p 7 X3,p 8 X 1 ,p 8 X 2 ,p 8 X3,p 9 Xi,p 9 X 2 ,p 9 X 3 ,piQX 1 , p 10 X 2 , p 10 X 3 } 

(4.3) 

has 1 ° 1 '! 1 (3_ 3 7)^ ! = ^0 elements. To find D (Subroutine |2|) we have to introduce a sequence 
of 30 dummy variables 

Y = {yi,y2,y3, 2/4, 2/5, 2/6, 2/7, 2/8, 2/9, 2/10, 2/11, yi2,yi3,yi4,yi5,yi6, 

2/17 , 2/18 , 2/19 , 2/20 , 2/21 , 2/22 , 2/23 , 2/24 , 2/25 , 2/26 , 2/27 , 2/28 , 2/29 , 2/30 } , (4.4) 
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expand the dot product G Y , and partition the terms in the resulting sum into equivalence 
classes. Two terms are equivalent if and only if they contain the same monomial in the 
variables {t\, t 2 , t 3 , £4, £5, ts, t 7 , x\, x 2 , x 3 }. Using the explicit form of the pi given in eq. 
( |4.2j ), we find 37 equivalence classes: 



{hxjyis, t 2 x\y X3 , t 3 x\y l3 , -t 5 xjy 16 , t 7 y l3 x\ + t 7 y 2% x\,t x x\x 2 yx^ t 2 x 1 x 2 yu ~ t 2 xix 2 y 4 , 
hxix 2 yi4, t 4 xix 2 y 22 - t 4 xix 2 y 25 , -t 5 xix 2 yn + t 5 xix 2 y 22 - t 5 xix 2 y 25 ,t G xix 2 y 22 - t 6 xix 2 y 25 , 
t 7 x x x 2 y 7 - t-jx^yxQ + t 7 xix 2 yu + t 7 x x x 2 y 22 - t 7 xix 2 y 25 + t 7 xix 2 y 29 , -t 2 x\y^ t 4 x\y 23 

-Uxly 2 6,t 5 xly 23 - t 5 xly 2G , t§x\y 23 - t§x\y 2 ^ t 7 y 8 xl - t 7 y n x 2 + hy 23 x\ - t 7 y 2 §x\, 

1 111 

-t\Xxx 3 y\Q + tiXixsym^^iXsyis - ~t 2 Xix 3 yi 3 xu, - ^t 3 x x x 3 yi + ^3x1x^4 

11 111 

--t 3 xix 3 y 7 + -t 3 xix 3 y 13 + t 3 x x x 3 yi 5 , -t A xix 3 yi 6 - -t 4 xix 3 yi 9 + -t 4 xix 3 y 22 

1 111 

--t4Xix 3 y 2 s,-t 5 xix 3 yis - -t 5 xix 3 y 22 xi4, + ■^t 5 x 1 x 3 y 25 xu,--^t & xix 3 y 25 ,t 7 x 1 x 3 y 15 

1 111 

+t 7 xix 3 y 30 , -hx 2 x 3 yn, -t 2 x 2 x 3 ye - -t 2 x 2 x 3 yuXi4, ~-^3X 2 x 3 y 2 + -t 3 x 2 x 3 y<r> 

11111 

--t 3 x 2 x 3 y s + -t 3 x 2 x 3 y 1A , -t 4 x 2 x 3 yi 7 - -t 4 x 2 x 3 y 20 + -t 4 x 2 x 3 y 23 + t A x 2 x 3 y 24l 

1 11 

-t 4 x 2 x 3 y 27 - -t4X 2 x 3 y 29 ,t 5 x 2 x 3 y 24 - t 5 x 2 x 3 y 27 - -t 5 x 2 x 3 y 23 xu + -t5X 2 x 3 y 2e xu, 

1 

t&x 2 x 3 y 24 - -t 6 x 2 x 3 y 26 - t 6 x 2 x 3 y 27 ,t 7 x 2 x 3 y g - t 7 x 2 x 3 y 12 + t 7 x 2 x 3 y 24 - t 7 x 2 x 3 y 27 , 

^hxjyi 2 , -^t 2 xjy l5 xi4, -^t 3 y 3 xj + -hy^xj - -t 3 y$xl + -t 3 yisx% -^m%x\ 

-~Uy2ixl + ~t4yux 3 - ^Uy 30 xl, ~t 5 xly 27 xi4 - ^S^IWxm, - j . (4.5) 

After setting the terms in each equivalence class to zero, we find that six of the resulting 
equations are obviously redundant. The remaining 31 equations, 

{0 = y 5 ,0 = yn,0 = yi2,0 = yi 3 ,0 = y 4 - y u ,0 = y 2 - y 5 + y$ - y u ,0 = y u , 

n 2/14X14 n on 1 n 2 ^ 15 

= — h 2/6, = yi - 2/4 + 2/7 - 2/13 - 22/15, = 2/3 - 2/6 + 2/9 - 2/15,0 = 2/13 , 

2 X14 

n n 1 o n n 2/22X14 2/25X14 „ „ 

= 2/15, = 2/10 + 22/15, = 2/16, = — - - h 2/18, = 2/22 - 2/25, = 2/25, 

= 2/17 - 2/22 + 2/25, = 2/23 - 2/26, = 2/8 ~ 2/11 + 2/23 ~ 2/26, = 2/24 ~ 2/27, 

2/26 2l/24 2j/27 

= 2/9 - 2/12 + 2/24 - 2/27, = 2/24 ^ 2/27, = 1 h 2/23 ~ 2/26, = 2/27, 

2 X14 X14 

= 2/16 - 2/19 + 2/22 - 2/28, = 2/13 + 2/28, = 2/17 - 2/20 + 2/23 + 22/24 - 22/27 - 2/29, 

= 2/7 - 2/10 + 2/14 + 2/22 - 2/25 + 2/29, = J/18 - 2/21 + 2/24 - 2/30, = 2/15 + 2/3o} , (4.6) 



constrain the j/j. After putting ( [4.6[ ) into matrix form, we can easily solve for the null space 
of the resulting matrix. We find that, as expected, the null space is one-dimensional: 

D = | (1,0, 0, 0, 0, 0, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1,0, 0, 0, 0, 0, 0, 0, 0, 1,0)}. (4.7) 
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We have now determined a basis for the scalar syzygies of P\: 

B 1 = D (4.8) 

To determine what syzygy of Po this element of So (Pi) corresponds to, we apply the same 
mapping that we used in Section ^ (formalized in Subroutine ||) in the example illustrating 
why computing a basis for a syzygy module with linear algebra is non-trivial. We partition 
the single element of D into ten distinct non-overlapping sequences of length three without 
disturbing the ordering of the vector that we started with 

{{1,0,0}, {0,0,0}, {-1,0,0}, {0,0,0}, {0,0,0}, {0,0,0}, {0,-1,0}, {0,0,0}, {0,0,0}, {0,1,0}} 

(4.9) 

(3) 

and then take the dot product of each sequence of three elements with Mj = {xi, X2,x 3 }: 

SyzList! = | (xi,0, -xi, 0, 0, 0, -x 2 , 0, 0, z 2 ) j • (4.10) 

It is easy to check that (xi, 0, — x%, 0, 0, 0, — X2, 0, 0, x 2 ) is indeed a syzygy of Po 

Now we increment d to two and begin our second pass through the While loop. This 
time, B\ ^ {0} and, therefore, the algorithm directs us to the upper branch of the If 
statement. M{ = {xi,x 2 ,x 3 }, P = {pi,p 2 ,P3,PA,P5,P6,P7,P8,P9,Pio}, and 



Pi = {PIXI , PlX 2 , P1X3, P2X1 , p 2 X 2 , P2X3 , P3X1 , P3X2, P3X3 , P4X1 , P4X2 ,P4X 3 ,p 5 X 1 ,p 5 X 2 ,P5X 3 , 
PqXi ,P6X 2 ,P6X3,P7Xl, P7X 2 , p 7 X 3 , p 8 X 1 , p 8 X 2 , p 8 X3, p 9 Xi, p 9 X 2 , PqX 3 , p 10 X 1 , p W X 2 , ^10^3 } 

(4.11) 

with ~^np~]iT^ = 30 elements. To determine P 2 , we need the set of monomials of degree 

two built out of {xi, X2, x 3 }, = {x\ , xix 2 , X1X3, x\, X2X3, x|}. As mentioned in Section 
^ for definiteness, we have chosen the lexicographical ordering for our sets of monomials. 
Let us say again that this choice is not by any means necessary. A different choice of 
ordering for the monomials would lead to a different representation of the same syzygy 
module; the number of linearly independent syzygies at each degree would be exactly the 
same. Continuing, we see that 

P2 = {pixl,PlXlX2,PlXlX3,pixl,piX 2 X3,Plxl,p 2 xl,P2XlX 2 ,p 2 XlX3,p 2 xl,p 2 X 2 X3,p 2 X3, 
P 3 xl,p 3 X 1 X2,p 3 X 1 X 3 ,p 3 xl,p 3 X 2 X 3 ,p 3 xl,p 4: xl,P4 : XlX2,P4X 1 X 3 ,P4xl,P4X2X 3 ,p i X 3 l , 
P5xl,P5XlX 2 ,p 5 XiX3,p 5 xl,P5X 2 X3,PBxl,Pexl,p 6 XiX 2 ,p 6 XiX3,Pexl,peX 2 X3,PQxl, (4.12) 
p 7 xl,p 7 X 1 X2,P7XlX 3 ,p 7 xl,p 7 X2X 3 ,p 7 xl,p 8 xl,p 8 X 1 X2,P8XlX 3 ,p 8 xl,psX2X 3 ,p 8 xl, 
p 9 xl,p 9 X 1 X2,P9XlX 3 ,p 9 xl,p 9 X2X 3 ,p 9 xl,p W xl,p W X 1 X2,PloXlX 3 ,pioxl,pi X2X 3 ,pioxl} 

has ^w|~]T! = 60 elements. Pi is the basis with respect to which the single element of 
B\ has components 

(1, 0, 0, 0, 0, 0, -1,0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1, 0, 0, 0, 0, 0, 0, 0, 0, 1,0) . (4.13) 
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The function a (Subroutine |l|) maps this vector to three new vectors in the larger vector 
space for which P2 is the standard basis. This map is carried out by extending the P\ by 
each of the three variables in turn. Applying a to 

(1, 0, 0, 0, 0, 0, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1,0, 0, 0, 0, 0, 0, 0, 0, 1, 0) 

= pixi -pzxi -p 7 x 2 +P10X2 , (4.14) 

we obtain 

CT {$\X\ - p 3 X X - p 7 X 2 + P10X 2 ) = {pix\ - p 3 x\ - p 7 X X X 2 + PlOX!X 2 , 

p\X\x 2 - p 3 x x x 2 - p 7 x\ + P\v,x\ y p\xxx 3 - p 3 x x x 3 - p 7 x 2 x 3 + P10X2X3} . (4.15) 

We can now read off the components of three scalar syzygies of P 2 which have their origin 
in the scalar syzygy of Pi that we computed on our first pass through the While loop. 
They are: 

A = |(1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
-1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0), (0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1, 0, 0, 
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1,0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
0, 1, 0, 0), (0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
0, 0, 0, 0, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1,0)}. (4.16) 

If we view the elements of A as the rows of a matrix, the result is already in row echelon 
form. It follows that C is simply A in matrix form. The first three columns of C are pivot 
columns, so G is P 2 with the first three entries replaced with zero: 

G = {d,Q,d : pixl,p l x 2 x 3 ,p l xl,p 2 x\,p 2 xix 2 ,p 2 x 1 x 3 ,p 2 xl,p 2 x 2 x 3 ,p 2 x 2 3 ,p 3 x\, 

P3XlX 2 ,p 3 XlX 3 ,p 3 xl, p 3 X 2 X 3 ,p 3 xl, P4X\ , p 4 Xi X 2 ,p±XlX 3 , p±x\ , P±X 2 X 3 ,Pax\ } p 5 x\ , 
P5XlX2,P5XlX3,p 5 xl,P5X2X3,p 5 xl,Pexl,PGX 1 X2,P6X 1 X 3 ,pexl,peX2X3,Pexl,p 7 Xi, (4.17) 
p 7 xi x 2 , p 7 xix 3 , p 7 x 2 2 , p 7 x 2 x 3 , p 7 xl , p 8 xl , p 8 X! x 2 , p 8 xix 3 , p 8 xl , p 8 x 2 x 3 ,p 8 xl, pgxj , 
p 9 x 1 x 2 ,p 9 xix 3 ,p 9 xl,p 9 x 2 x 3 ,p 9 xl,p 10 xl,p w xix 2 ,p w xix 3 ,pioxl,p w x 2 x 3 ,pi xl} . 

Making this replacement will prevent Subroutine || from rediscovering syzygies that have 
their origin in the syzygy of Pi that we found on our first pass through the While loop. On 
this pass through the loop, we refrain from stepping into Subroutines Q and ^ since their 
functionality should already be quite clear from our first pass through. Applying SSyz to 
G, we find 

D = SSyz(G) = I (0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
0, 0, -1, ^,0, -4, -1,0, -2, ^,0, -2, -1, 0, 0, 0, 0, -2, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0) , (0, 0, 0, 0, 
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1 Yi 4 1 Y14 111 1 1 Y14 

--,^,0,0, --,0,0,^,0, 0,0, 0,0, -,0,0, 0,0, 0,0, 0,0, --,0,0, 0,0,0,-,-^, 
0,2,1, -1,1, -^,0,1,^,0,0, 0,0, 1,0, 0,0, 0,0, 0,0, 1,0,0) I (4.i 8 ) 
for the new scalar syzygies of P 2 . Since D ^ {0}, 

B 2 = C U D = |(1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
0, 0, 0, 0, -1,0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0), (0, 1,0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
-1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1,0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
0, 0, 0, 0, 1, 0, 0), (0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 

0, 0, 0, 0, 0, 0, 0, -1,0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0), (0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 



0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1, £1,0, -4, -1, 0, -2, £1,0, -2, -1,0, 0, 0, 

0,-2, 0,0, 0,0, 0,0, 0,0, 1,0), (0,0, 0,0, -i,^,o, o,-i,o, 0,^,0, -1,-1, 0,0, 0,0, 1,0, 

0,0,0,0,0,0,0,-1 0,0,0,0,0,1 -£1,0,2,1 -1,1 -*H 0,1,1 0,0,0,0,1,0,0,0,0,0,0, 

0,1,0,0)}. (4.19) 

The two new scalar syzygies of P 2 discovered by SSyz map, via ToSyz, to two new syzygies 
of Po of degree two, linearly independent of (x±, 0, — x\, 0, 0, 0, — x 2 , 0, 0, x 2 ): 

SyzList 2 = ToSyz(D) = j (0, 0, 0, 0, 0, I23X14 - x 2 x 3 , ^xjxu - xix 3 - 2x 2 x 3 

-4xix 2 , -2xix 2 - xix 3 , -2xix 2 , x 2 x 3 ) , (l^Xu - ^p, 1^3X14 - -^p, 

1 X1X3 X\X 2 1 X 2 X 3 1 2 1 2 2 

~-jXix 2 — , — — , --x 2 x 3 , — ^3X14, -^£3X14 ~x 2 + 2xix 2 

x 3 x 2 x x x 3 Xl x 3 2 \ \ lA on N 

H — y~ + ~^2~ ,XlX2 ~Y~ ' XlX2 ' x 2j } ■ ( 4 - 20 ) 

It is again simple to check that the elements of SyzList 2 are in fact syzygies of Po- If we 
did not have a priori knowledge of the syzygy module, we would have to continue on and 
pass through the While loop several more times (as explained in |l0| , A = 6 for the planar 
massless double box) to search for more potentially useful linearly independent syzygies. 15 
Since we do not expect to find further solutions, we can exit the While loop after just two 
iterations and collect the results: 

2 f 11 

(J SyzList d = < (xi, 0, -x 1 ,0, 0, 0, -x 2 , 0, 0, x 2 ) , (o, 0, 0, 0, 0, -X3X14 - x 2 x 3 , -xjxu 

d=i 

-X1X3 - 2x 2 x 3 - 4xix 2 , -2xix 2 - xix 3 , -2x x x 2 , x 2 x 3 ) , (l^xu - -^p, -xlxu 



15 It is important to realize that, in general, D — {0} for d — i does not imply that D = {0} for d — i + 1, 
though this may turn out to be true in practice for physically motivated examples. 
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XlX 3 1 XiX 2 1 X 2 X 3 1 2 1 2 2 

£-» -2 X i x 2 -g-, -^2*3, — ^3X14, -^^3X14 - £ 2 + 2si*2 

. x 3 x 2 xix 3 xix 3 2 \ 1 

H — jj— + ~^-,xix 2 + ~^-,xix 2 ,x 2 J S . (4.21) 

Before leaving this section, we translate the above result back into the usual language used 
to describe the planar massless double box: 

2 , 

|J SyzList d = m ■ & 4 , 0, -4 • fc 4 , 0, 0, 0, • A*, 0, 0, 4 ■ fci) , (o, 0, 0, 0, 0, -s? 2 Xi4 (4.22) 
1 2 

-«2 • &i «i2, 2 S i2Xi4 -£\-k4 S12 - 24 • A;i S12 - 44 • k^h • h, -24 ■ fc 4 4 • fci - 4 • fc 4 si 2 , 

-24 • fc 4 4> • fci, 4 ' Kl«12 J, (^12X14 ^ ' g S 12Xl4 1 ' ~2 ' ' 

4-&4S12 £l-k 4 £ 2 -ki 1 4-^lSl2 1 2 1 2 C/> j \2 

4 , 2 > -4^2 • «i S12, 2 4*12X14, -^s 12 xi4 - (4 ■ ki) 

+24 • fc 4 4 • fci + + ^p,4 ■ M2 ■ fci + ■ M2 • fci, (4 • fci) 1 

5. Conclusions 

This paper continued the program of research initiated in ]l(J by Gluza, Kajda, and 
Kosower (GKK). In reference jl(J], GKK proposed a novel reduction scheme for multi- 
loop integrals guaranteed to produce unitarity-compatible integral bases free of doubled 
propagator denominators. In this work we took a fresh look at the computationally in- 
tensive part of their procedure, the generation of complete sets of unitarity-compatible 
integration by parts relations. Drawing upon some of the ideas in reference we found 
an attractive alternative to the GKK approach which completely avoids the use of Grobner 
bases. In fact, we showed in Section ^ that our solution, Algorithm [l], can be described in 
terms of simple linear algebra. 

One shortcoming of the present paper is that we cannot yet claim to have fully opti- 
mized Algorithm |l|. Even if we excise the trivial redundancies introduced into our pseudo- 
code for the sake of clarity, there are many features of the algorithm which might benefit 
from further optimization. For example, as should already be clear from the non-trivial 
example discussed in Section ||, the matrices produced by Algorithm p] are typically quite 
sparse and, so far, we have made no attempt to exploit this feature of the problem. Fur- 
thermore, it seems likely that a high-level implementation in Mathematica (such as the 
one written by the author) will ultimately prove insufficient for research purposes. Besides 
the fact that, as a general rule, Mathematica runs very slowly, it is not at all clear that 
Mathematica exploits the best available algorithms for row reducing matrices; if our pre- 
liminary experimentations are any guide, it appears that Mathematica manages system 
resources rather poorly. Fortunately, we anticipate that the brevity and simplicity of Al- 
gorithm [l] will make it possible to optimize and then effectively implement at a lower level 
in C++ or Fortran. We are excited by this prospect and hope to pursue a project along 
these lines in the near future. 
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